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ELASTIC-PLASTIC  ANALYSIS  OF  CLOSED-END  TUBES 


ABSTRACT 

approximate  solution  is  derived  for  the 
stresses  in  a  circular,  closed-end,  unrestrained 
tube  subjected  to  an  internal  pressure  of  suffi¬ 
cient  magnitude  to  Induce  plastic  flow.  The 
Mises  yield  criterion  is  assumed  in  the  plastic 
region  of  the  tube.  The  radial  and  circumfer¬ 
ential  stresses  are  found  to  be  in  good  agreement 
with  those  for  the  plane  stress  case.  Results 
are  shown  graphically. 
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CONCLUSIONS 


By  assuming  a  logarithmic  distribution  for  the  longitudinal  stress 
in  the  plastic  portion  of  the  tube  wall  and  comparing  these  results  with 
other  numerically  obtained  solutions  for  the  pressurized,  closed-end 
tube,  it  has  been  found  that  very  little  difference  in  the  radial  and 
tangential  stress  distribution  is  observed.  Of  course,  the  end- force 
must  be  equilibrated  by  the  longitudinal  stresses  acting  across  the 
tube  thickness,  but  as  long  as  this  condition  is  met  it  appears  that 
any  assumed  distribution  for  the  O'*  stress  in  the  plastic  region,  has 
little  effect  upon  the  aforementioned  radial  and  tangential  stresses. 

One  point  which  should  be  mentioned  is  that  Che  Mises  yield 
criterion  was  considered  to  govern  the  behavior  of  the  tube  material  in 
the  plastic  state.  Hence,  the  fact  that  the  CT*  and  CT#  stresses 
do  not  differ  greatly,  regardless  of  what  end  conditions  for  the  tube 
are  considered,  is  not  totally  unexpected.  Particular  reference  is 
made  to  the  solution  of  this  same  problem  wherein  the  Tresca  yield 
criterion  is  utilized.  In  this  case,  it  is  found  that  the  end  condition 
has  no  Influence  upon  (J^  and  0?  • 

Comparison  of  these  same  stresses  with  tho  plane  strain  solution1 
indicates  approximate  agreement.  This  is  not  surprising  in  view  of  the 
fact  that  it  oan  be  readily  shown  that  the  magnitude  of  the  end  restraint, 
which  is  required  to  maintain  the  plane  strain  condition,  approaches  the 
end  force  produced  in  the  closed-end  tube  for  the  case  of  complete 
plasticity*. 

'i 

A  fair  agreement  with  the  plane  stress  results  Is  also  obtained  . 

In  fact,  for  the  plane  strain,  plane  stress  and  closed-end  conditions, 
it  is  interesting  to  note  that  very  little  difference  is  observed  in  the 
distribution  and  magnitude  of  the  radial  and  tangential  stresses  for  the 
case  of  initial  yielding  at  the  bore  surface. 

In  comparing  the  magnitude  of  the  longitudinal  stress  CJi  ,  as 
developed  herein,  with  the  values  of  0#  obtained  by  use  of  the 
Prandtl-Reuss  stress  strain  relations  for  a  perfectly  plastic  material, 
large  discrepancies  are  observed.  Furthermore,  it  is  generally  agreed 
that  use  of  the  Prandtl-Reuss  relation  leads  to  more  accurate  stress 
predictions;  however,  the  mathematical  difficulties  encountered  are 
more  formidable. 
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CONCLUSIONS  (CONT.) 

In  any  event,  useful  stress  predictions  for  0^  and  C 1?  can  be 

obtained  by  the  methods  which  are  outlined  in  this  report. 
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INTRODUCTION 


The  problem  is  one  of  a  thick-wal led ,  circular,  closed  end, 
unrestrained  tube  subjected  to  an  internal  pressure  of  sufficient 
magnitude  to  induce  plastic  flow.  The  plastic  deformation  occurs 
first  at  the  inner  surface  of  the  tube  and,  as  the  pressure  is 
increased,  the  plastic  region  continues  to  grow  until  the  entire 
tube  becomes  plastic.  A  state  of  contained  plasticity  is  assumed. 
The  tube  material  is  considered  to  be  homogeneous,  Isotropic,  com¬ 
pressible,  and  perfectly  plastic. 

To  describe  the  pressure  induced  stresses  and  strains  in  the 
tube,  a  system  of  cylindrical  coordinates  (  ,  6 , 2.  )  is  employed, 

where  the  2  -axis  coincides  with  the  axis  of  the  tube.  All  stresses 
and  strains  are  assumed  to  be  independent  of  0  and  2  . 


FORMULATION 

The  principal  directions  of  stress  at  a  generic  point  of  the 
tube  are  radial,  circumferential  and  axial,  and  the  only  equation  of 
equilibrium  which  remains  to  be  satisfied  is^: 


^  ^  _L  Or  Qo  -  q  (1) 

dr  r 


We  denote  the  inner  surface  of  the  tube  by  I f  <=  &.  ,  the  outer 
surface  by  If  =  b  ,  and  the  interface  between  the  elastic  and  plastic 
region  by  p  “  /  •  We  denote  by  Or*,  G#*  ,  and  0^p  ,  0JP 

the  stresses  in  the  elastic  and  plastic  regions,  respectively.  The 
Mises  yield  criterion  is  assumed  and  hence,  throughout  the  plastic 
region  A  ^  f  £  f  t  A  <  j5  $  b  »  the  following  condition  must  be 
satisfied: 


(Vrr-  07'/  +  (  07'  -  <J//+  ( <37"  -  07' /=  6  (2) 


where  k 


is  the  yield  stress  in  simple  shear  for  the  material. 
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The  solution  for  the  elastic  region  P$  p^;  \q  ,  A.  $  £  $  lo  ,  la 
given  by  the  Lam&  equations: 


c;e=  A  +  sljuV  -  2fJL6/r' 
o;c=  A  +  +  «2/u  3/r* 

o;*=  A  <3> 


The  boundary  conditions  to  be  satisfied  are: 


i.  07=  -  P  ;  r  -  <x 

Z.  07  *  0  j  r  “  10 


3. 

4. 


o;e-  07  ,  cT~  a/  ,  0/ -  07p  , 
=^Tr[c7ir^=27rf  WrAir  + 


r  =  p 

■  b 

C V*xLr  > 

V 


The  last  condition  states  that  the  axial  stress  integrated  over  the 
thickness  of  the  tube  equilibrates  the  force  acting  on  the  closed  end 
of  the  tube.  These  conditions  are  sufficient  to  determine  the  constants 
of  integration  and  establish  a  relationship  between  p  and  . 

To  determine  the  value  of  p*'  P  for  incipient  plastic  flow, 
i.e.  ?=■&.)  we  apply  the  boundary  conditions  1,  2,  and  4  to 
eauations  (3)  and  find  the  constants  A,  B,  D  in  terms  of  p 
This  yields; 


A- j££- 
I*'-  <C 


\5  * 


> 


p*fC  lo 


D  -  o  . 
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r  ~  a.  , 


We  next  apply  the  Mises  yield  condition,  equation  (2),  setting 
and  soLve  for  p*  ; 


*=  tt  ( bv-  Q 


To  determine  the  plastic  stress  state  for  P  >  A.  ,  we  first 
solve  (2)  for  :  * 

a*P  =  +  ^\)  +Sr-To <4> 

The  ambiguous  sign  (  ±  )  in  (4)  is  determined  by  using  the  results 
for  incipient  plastic  flow,  evaluating  the  stresses  for  J"  sa.  ft. ,  and 
substituting  for  the  stresses  in  (4).  It  follows  that  the  plus  sign 
should  be  used  in  (4). 

We  introduce  the  notations: 


-  37.  « 

afe  j  — 

2k 

(ft  = *  -A- 

D 

; 

=  JEL  p  . 

2k  j  V 

JL 

2k 

The  equation  of  equilibrium  and 

equation  (4)  become: 

Substituting  for  0"^  P  from  (4)  into  (1),  and  letting 
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/U  07  ?  we  obtaln 


M 

+  -M 

zr  1 

r  - 

-  ij  3(  \  -  ^l)  j 

(5) 

Letting 

M- 

^3d-<£)  , 

(6) 

ar 

A  3T 

+  JSL  -  0 

(7) 

d  r 

ar 

two  equations  arc  obtained  which  arc  immediately  integrable,  Integrat¬ 
ing  (b),  wc  find; 


1 


(8) 


where  C|  is  a  constant  ol  integration.  Substituting  for  ^  from 
(8)  into  (7)  and  integrating,  we  obtain:  11 


<37 


P 


(9) 


where  is  a  constant  of  integration.  From  (8)  and  (A),  we  now  find: 


+ 


+•  C*.  (10) 
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and 


These  six  equations  involve  the  six  unknown  quantities  ft  >  ®  > 

C,  t  Cj.  and  .  The  first  four  of  these  equations  are  used 

to  eliminate  D  and  C*.  from  the  last  two.  After  simplification 

and  the  introduction  of  the  notations 
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these  two  equations  take  the  form: 

smC,  l^'E\3b'-7)]smf  +  48  ^ 

+  3  £(3^-11)  J1" -b  S'?  fcT’C  E*’-' ()]  +  aagb  p 

+dos?,[3[(3^-«)jVs7Hvc^-0]»«7  +  6finr  ^ 

+  fy[»(rt-is^)r+i'^stf''',3“*?}  “  0  tl8) 

sm?j"-P(3?-7W*w?-3(»^*Vpv1)e»»f  +  ia^j  4-  la^p 
+  toi  ?,  [-  3  (3  *">  f )  sm  |  +•  (5  (3  f  -  7  s')  «>  ?  ^ 1 ***  £  }  -  °(19) 


Equations  (18)  and  (19)  are  then  solved  for_  Slfi  C,  and 
COS^i  i n  terms  of  p  and  f  (  p  and  f  ),  and  the 
relationship  between  p  and  j3  is  then  determined  from  the  identity 

cos'-c,  +  Sin'-Q,  =  \ 


The  resulting  expression  is: 


Jl 

(20)J 


p  a Ei*^j |aVsi  E%[-8if3(fv-  T}’)s\nJ-3Z ~B\osJ + 3 J - (>?  1\  llfcVte] 

where 

74bl^  =  ^3(Et+0f‘l+  nt^SW?-l6lI5^-0?'+a 

+  f3  [3f f-  6Jt'bV+  19  b*'  +  /6  b^jj 


Ci>S 

(21) 


If  we  let  then  ?-!>  .P  =  ^  fi<v\  J  =»•  o  and  (20) 

reduces  to  the  previously  determined  result  for  incipient  plastic  flow: 

P*=  ^--1 
'  a  ¥ 
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For  ,  the  entire  tube  is  in  the  plastic  state  and  the  pressure 

corresponding  to  this,  "n  a.  p**  ,  is  found  from  (20)  by  letting 

p  *  *  _  X  fjT Oft  —  eas'£)+Mi  -  l  )  si»\'b} 


57  fr*  |  _  r-i. 

j  o  +-  I  —  2  b  cos  b  v 


The  values  of  the  various  constants  are  given  by: 

C05  c,  =  y[7f^]  |*^T(|  sp  -  Bl;  smp  +  3 ( pV  wsp  -  6 o| 

cl=[7^][4(?f+,3Ca)sinp-/tf^f-2B>5p-ir5C3?t6fBV,^T5¥jj 

A  |ia^ ?V+  ^)sinp-4/3(4f  -  3  &)usj-  ftfCs?-  7  Wj 


B 


—  a.  j 

_  b  1 

r  p  “1 

^  i 

1.76  B%AJ 

jiapVnf -it&f’usJ-fsjYsf-  lit';  j 


Introducing  the  notation  r  *•  %.  >  IT  *  ^  Cvv  T  we  find 
the  stresses  for  the  plastic  region  y.  ^  ^ 
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fcf« 


74E"A 


j"4^  (*f%-  nS'jsmCF-  f  )4-  latff  j'"£Ubos,t'r-f) 
-/^E  (WsiVr  -i'  \T5c.osr)+  4(s  f\  13 k1)  5'm^ 

-  / 6  f3  [f^-S^ccsf  -  ft  (jJ*.  6  f  E\  H  bfjj. 

£■*  f- 11?  E>in  ( r- J)+aif3 

-8  Ev(^ ***7 +6|f3aM'r.)  +  4  (3 |V  13 

-  /  6  lfT  ( f-  a  5V)  cmJ  -  tf5  (3  Jl  if  &'+• 1  <1  E+j] 

^ P=  [wS**  {'  ^ J?  - '’  5>  '*< ?-  ?>+ <*  ®tf*+  n  e)»s(  7- 1 ; 
- 4El(sm7  +  s-if3 tos  r)  -h4C3f  +  i3¥^ 

-/4  (T(  f  -  a  E‘)  cos?  -  fj (3  4  f  B'+  II  b’jj 

given1  by ;  tl'C  ClaStic  re8ion  f  *  r  £  fc>  »  A  £  f  £  fc>  ,  the  stresses 


are 


of*  [-!_- 

r  1_  7i  FA 


lasing  -  IC^3  ccsJ^\[l(3fx-i^hx)Jfxilz^J 


°*  *  =  [7^3]  a  sin  5  -  /fepeos  J  -  J3  (3  f- 11  f)jfCr%p 
°**  ”  [7^]  [ia(f'++BVsinf  -  4  (5  (4f--  3  5*)  cosj 

-  If?  fX  3 j>*_  7  B'jJ 
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RESULTS 


Numerical  results  were  obtained  for  the  case  b  =*  v«.  53  2. 

1.0  ,  USL,  /.*,  /.6  J  1.2,  2.6 

The  equations  for  p  ^  p/3^  and  the  various  stresses  were  coded 

for  the  IBM  650  digital  calculator.  The  results  are  shown  graphically 
in  figures  one  through  five. 
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